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Dynamical  Theory  for  Treating  the  Motion  of 
Cold  and  Warm  Fronts  in  the  Atmosphere 

by  J.  J,  Stoker 


Sec.  1.   Introduction 

One  of  the  most  characteristic  features  of  the  motion 
of  the  atmosphere  In  middle  latitudes  and  also  one  which  is 
of  basic  importance  in  determining  the  weather  there  is  the 
motion  of  wavelike  disturbances  which  propagate  on  a  discon- 
tinuity surface  between  a  thin  wedge-shaped  layer  of  cold 
air  on  the  ground  and  an  overlying  layer  of  warmer  air.   In 
addition  to  a  temperature  discontinuity  there  is  also  in 
general  a  discontinuity  in  the  tangential  component  of  the 
wind  velocity  in  the  two  layers.   The  study  of  such  phenom- 
ena was  initiated  long  ago  by  Bjerknes  and  Solberg  [2] 
and  has  been  continued  since  by  many  others.   In  considering 
wave  motions  on  discontinuity  surfaces  it  was  natural  to  be- 
gin by  considering  motions  which  depart  so  little  from  some 
constant  steady  motion  (in  which  the  discontinuity  surface 
remains  fixed  in  space)  that  linearizations  can  be  per- 
formed, thus  bringing  the  problems  into  the  realm  of  the 
classical  linear  mathematical  physics.   Such  studies  have 
led  to  valuable  insights,  particularly  with  respect  to  the 
question  of  stability  of  wave  motions  in  relation  to  the 
wave  length  of  the  perturbations.   (The  problems  being  lin- 
ear, the  motions  in  general  can  be  built  up  as  a  combina- 
tion, roughly  speaking,  of  simple  sine  or  cosine  waves  and 
it  is  the  wave  length  of  such  components  that  is  meant  here, 
cf.  Haurwitz  [5],  p.  2-34.)  One  conjectiore  is  that  the  cy- 
clones of  the  middle  latitudes--probably  the  most  striking 
single  meteorological  phenomenon  there--are  initiated  because 

Numbers  in  square  brackets  refer  to  the  bibliography  at 
the  end  of  this  report. 
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of  the  occurrence  of  such  imstable  waves  on  a  discontinuity 
surface. 

A  glance  at  a  weather  map,  or,  still  better,  an  exami- 
nation of  weather  maps  over  a  period  of  a  few  days,  shows 
clearly  that  the  wave  motions  on  the  discontinuity  surfaces 
(which  manifest  themselves  as  the  so-called  fronts  on  the 
ground)  develop  amplitudes  so  rapidly  and  of  such  a  magni- 
tude that  a  description  of  the  wave  motions  over  a  period 
of,  say,  a  day  by  a  linearization  seems  not  feasible  with 
any  accuracy.   The  object  of  the  present  report  is  to  make  a 
first  step  in  the  direction  of  a  nonlinear  theory,  based  on 
the  exact  hydrodynamical  equations,  for  the  description  of 
these  motions  that  can  be  attacked  by  numerical  or  other 
methods.   The  reason  that  no  such  treatment  seems  to  have 
been  made  up  to  now  is  not  that  meteorologists  have  failed 
to  see  the  desirability  of  doing  so,  but  rather  that  the 
difficulties  to  be  overcome  in  obtaining  even  numerical  so- 
lutions are  really  formidable.   No  claim  is  made  that  the 
problem  is  solved  here  in  any  general  sense.   What  is  done 
is  to  start  with  the  general  hydrodynamical  equations  and 
make  a  series  of  assumptions  regarding  the  flow;  in  this  way 
a  sequence  of  three  nonlinear  problems  (we  call  them  Prob- 
lems I,  II,  III),  each  one  furnishing  a  consistent  and  com- 
plete mathematical  problem,  is  formulated.   In  this  way  one 
can  see  the  effect  of  each  additional  assumption  in  simpli- 
fying the  mathematical  problem.   The  first  two  problems  re- 
sult from  a  series  of  assumptions  which  would  probably  be 
generally  accepted  by  meteorologists  as  reasonable,  but  un- 
fortunately even  Problem  II  is  still  unmanageable  in  terms 
of  numerical  analysis.   Further,  and  more  drastic,  assump- 
tions lead  to  a  still  simpler  Problem  III  which  is  formu- 
lated in  terms  of  three  first  order  partial  differential 
equations  in  three  dependent  and  three  independent  variables 
(as  contrasted  with  eight  differential  equations  in  four  in- 
dependent variables  in  Problem  I).   The  three  differential 
equations  of  Problem  III  are  probably  capable  of  yielding 
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reasonably  accurate  approximations  to  the  frontal  motions 
under  consideration,  but  they  are  still  too  difficult  to 
deal  with,  even  numerically,  principally  because  they  in- 
volve three  independent  variables  ^'"'- -such  equations  are 
well  known  to  be  beyond  the  scope  of  even  the  most  modern 
digital  computing  machines.   Consequently,  still  further 
simplifying  assumptions  must  be  made  in  order  to  obtain  a 
theory  capable  of  yielding  some  concrete  results  through 
calculation. 

At  this  point,  two  different  approaches  to  the  problem 
have  been  devised.   One  of  them,  by  \¥hltham  (and  which  is 
discussed  by  him  in  this  report),  deals  rather  directly  with 
the  three  differential  equations  of  Problem  III.   Two  of 
these  equations  are  essentially  the  same  as  the  classical 
gas  dynamics  equations  In  two  Independent  and  two  dependent 
variables  if  one  regards  the  third  Independent  variable  as  a 
parameter  and  makes  one  additional  and  rather  reasonable  as- 
siimptlon.   These  two  equations --which  refer  to  motions  in 
vertical  planes --can  therefore  be  Integrated.   Afterwards 
the  transverse  component  of  the  velocity  is  found  by  inte- 
grating a  linear  first  order  partial  differential  equation. 
In  this  way  a  quite  reasonable  qualitative  description  of 
the  dynamics  of  frontal  motions  can  be  achieved  which  is  in 
good  agreement  with  many  of  the  observed  phenomena.   How- 
ever, this  theory  has  a  disadvantage  in  that  it  does  not 

The  work  of  Freeman  [3,4]  is  based  on  a  theory  which 
could  be  considered  as  a  special  case  of  Problem  III  in 
which  the  Coriolis  terms  are  neglected  and  the  motion  is 
assumed  at  the  outset  to  depend  on  only  one  space  varia- 
ble and  the  time.   The  idea  of  deriving  the  theory  re- 
sulting in  Problem  III  occurred  to  the  author  while  read- 
ing Freeman's  paper,  and,  indeed.  Freeman  indicates  the 
possibility  of  doing  so.   The  guiding  idea  is  to  make  use 
oi  the  analogue  of  the  long-v/ave  shallow  water  theory  for 
gravity  waves,  and  this  idea  seems  to  have  occurred  to  a 
number  of  meteorologists  (in  particular,  see  Abdullah 
Llj)._  However,  the  present  paper  seems  to  be  the  first 
in  which  essentially  three-dimensional  motions  are  stud- 
ied. 
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permit  a  complete  niamerical  integration  because  of  a  pecu- 
liar difficulty  at  cold  fronts.   (The  difficulty  stems  from 
the  fact  that  a  cold  front  corresponds  in  this  theory  to 
what  amounts  to  the  propagation  of  a  shock  into  a  vacuum--a 
mathematical  impossibility.   If  one  had  a  means  of  taking 
care  of  turbulence  and  friction  at  the  ground,  it  would  per- 
haps be  possible  to  overcome  this  difficulty.)   Neverthe- 
less, the  good  qualitative  agreement  with  the  observed  phe- 
nomena is  an  indication  that  the  three  differential  equa- 
tions fiornishing  the  basic  approximate  theory  from  which  we 
start--i.e.  those  of  our  Problem  III--have  in  them  the  pos- 
sibility of  furnishing  reasonable  results. 

The  author's  method  of  treating  the  three  basic  differ- 
ential equations  is  quite  different  from  that  of  VJhitham, 
but  it  \inf ortunately  involves  a  further  assumption  which  has 
the  effect  of  limiting  the  applicability  of  the  theory.   The 
guiding  principle  was  that  differential  equations  in  only 
two  independent  variables  should  be  obtained  in  order  to 
make  niAmerical  computations  feasible.   On  the  other  hand, 
the  nuaber  of  dependent  variables  need  not  be  so  ruthlessly 
limited.   Finally,  it  is  highly  desirable  to  obtain  differ- 
ential equations  of  hyperbolic  type  in  order  that  the  method 
of  characteristics  can  be  employed  as  an  aid  in  computing 
solutions.   These  objectives  can  be  attained  by  making  quite 
a  few  further  simplifying  assumptions  with  respect  to  the 
mechanics  of  the  situation.   The  result  is  what  might  be 
called  Problem  IV.   What  remains  of  the  original  hydrod3niam- 
ical  equations  is  still  complicated  enough,  however,  to  make 
one  feel  that  if  such  a  simplified  theory  fails  to  mirror 
the  observed  facts  at  least  roughly,  then  the  problem  of 
calculating  the  unsteady  motion  of  cold  and  warm  fronts  will 
present  a  formidable  challenge  to  numerical  analysis  in  the 
general  field  of  nonlinear  partial  differential  equations. 
The  theory  formulated  in  Problem  IV  is  embodied  in  a  system 
of  four  nonlinear  first  order  partial  differential  equations 
of  hyperbolic  type  in  four  dependent  and  two  independent 
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variables.   A  numerical  integration  of  these  equations  is 
being  carried  out  (by  finite  differences  applied  to  the 
equations  in  characteristic  form),  but  the  labor  of  inte- 
grating the  equations  is  so  great  that  only  meagre  results 
are  so  far  available. 

In  section  2  the  sequence  of  approximate  theories  for- 
mulated as  Problems  I,  II,  III  will  be  derived  from  the  hy- 
drodynamical  equations  by  adding  successively  more  and  more 
hypotheses  of  a  physical  character.   In  section  3  the  au- 
thor's method  of  treating  Problem  III--whlch  should  be  re- 
garded as  the  basic  problem--approximately  will  be  derived. 
This  yields  Problem  IV.   In  section  4  the  characteristic 
form  of  the  differential  equations  of  Problem  IV  will  be  ob- 
tained and  conclusions  drawn  from  them.   This  is  foliov/ed  by 
the  report  of  G.  \Vhitham  giving  his  method  of  treatment  of 
the  differential  equations  of  the  basic  Problem  III.   The 
theory  of  G.  Whitham  is  in  some  Virays  more  satisfactory  than 
that  of  the  author  since  it  yields  important  qualitative  re- 
sults of  the  kind  actually  observed,  and  does  that,  in  addi- 
tion, without  the  necessity  to  make  numerical  calculations. 
However,  it  is  unsuitable  for  describing  the  early  stages  of 
frontal  motions,  and  this  suggests  that  the  method  of  the 
author  might  be  used  for  the  initial  stages  followed  by  the 
method  of  Whitham  for  the  later  stages  of  the  motions. 

Sec.  2.   Derivation  of  the  basic  approximate  theory 

In  the  present  section  the  details  of  the  derivation  of 
the  basic  approximate  theory  will  be  given,  but  it  seems 
worth  while  to  sketch  out  in  general  terms  the  underlying 
ideas  and  motivations  first  by  way  of  orientation. 

To  begin  with,  a  certain  steady  motion  (called  a  sta- 
tionary front)  is  taken  as  an  initial  state,  and  this  con- 
sists of  a  iiniform  flow  of  two  superimposed  layers  of  cold 
and  warm  air,  as  indicated  in  Figure  1.   The  z-axis  is  taken 
positive  upward  and  the  x,y-plane  is  a  tangent  plane  to  the 
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earth.   The  rotation  of  the  earth  is  to  be  taken  into  account 
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but,  for  the  sake  of  simplicity,  not  its  sphericity--a  com- 
mon practice  in  dynamic  meteorology.   The  coordinate  system 
is  assumed  to  be  rotating  about  the  z-axis  with  a  constant 
angular  velocity  CI   =  60  sin  ^,  with  co  the  angular  velocity 
of  the  earth  and  50  the  latitude  of  the  origin  of  our  coor- 
dinate system.   (The  motivation  for  this  is  that  the  main 
effects  one  carr-s  about  are  found  so  long  as  the  Coriolis 
term  is  included,  and  that  neglect  of  the  curvature  of  the 
earth  has  no  serious  qualitative  effect.)   As  indicated  in 
Figure  1,  the  cold  air  lies  in  a  wedge  under  the  warm  air 
and  the  discontinuity  surface  between  the  two  layers  is  in- 
clined at  angle  OC  to  the  horizontal.   The  term  "front"  is 
always  applied  to  the  intersection  of  the  discontinuity  sur- 
face with  the  grovwid,  and  in  the  present  case  we  have  there- 
fore as  Initial  state  a  stationary  front  running  along  the 
X-axis.   The  wind  velocity  in  the  two  layers  is  parallel  to 
the  X-axis  (otherwise  the  discontinuity  surface  could  not  be 
stationary),  but  it  will  in  general  be  different  in  magnitude 
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and  perhaps  even  opposite  in  direction  in  the  two  layers. 
The  situation  shown  in  Figure  1  is  not  uncommon.   For  in- 
stance, the  X-axis  might  be  in  the  eastward  direction,  the 
y-axis  in  the  northward  direction  and  the  warm  air  would  be 
moving  in  the  direction  of  the  prevailing  westerlies.   The 
origin  of  the  cold  air  at  the  ground  is,  of  course,  the 
existence  of  the  cold  polar  regions.   We  shall  see  later 
that  such  configurations  are  dynamically  correct  and  that 
the  angle  QC  is  uniquely  determined  (and  quite  small,  of  the 
order  of  -|°)  once  the  state  of  the  warm  air  and  cold  air  is 
given.   (The  discontinuity  surface  is  not  horizontal  because 
of  the  Coriolls  force  arising  from  the  rotation  of  the 
earth. ) 

We  proceed  next  to  describe  what  is  observed  to  happen 
in  many  cases  once  such  a  stationary  front  starts  moving. 
In  Figure  2  (see  next  page)  a  sequence  of  diagrammatic 
sketches  is  given  which  indicate  in  a  general  way  what  can 
happen.   All  of  the  sketches  show  the  intersection  of  the 
moving  discontinuity  surface  (cf.  Figure  l)  with  the  ground 
(the  x,y-plane  with  the  y-axis  taken  northward,  the  x-axls 
taken  eastward).   The  shaded  area  indicates  the  region  on 
the  ground  covered  by  cold  air,  while  the  unshaded  region  is 
covered  at  the  ground  by  warm  air.   Of  course,  the  cold  air 
always  lies  in  a  thin  wedge  under  a  thick  layer  of  warm  air. 
In  Figure  2a  the  development  of  a  bulge  in  the  stationary 
front  toward  the  north  is  indicated.     Such  a  bulge  then 

What  agency  serves  to  initiate  and  to  maintain  such  mo- 
tions appears  to  be  a  mysterj''.   Naturally  such  an  impor- 
tant matter  has  been  the  subject  of  a  great  deal  of  dis- 
cussion and  speculation,  but  there  seems  to  be  no   consis- 
tent view  about  It  among  meteorologists.   In  applying  the 
theory  derived  here  no  attempt  is  made  to  settle  this 
question:   we  simply  take  our  dynamical  model,  assume  an 
initial  condition  which  in  effect  states  that  a  bulge  of 
the  kind  just  described  is  initiated,  and  then  study  the 
subsequent  motion  by  Integrating  the  differential  equa- 
tions subject  to  the  appropriate  initial  conditions.   How- 
ever, if  the  approximate  theory  is  really  valid,  such 
studies  might  perhaps  be  used,  or  could  be  modified,  in 
such  a  way  as  to  throw  some  light  on  this  Important  and 
vexing  question. 
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Stages  in  the  motion  of  a  frontal  disturbance 
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deepens  and  propagates  eastward  with  a  velocity  of  the  order 
of  500  miles  per  day.   It  now  becomes  possible  to  define  the 
terms  cold  front  and  warm  front .   As  indicated  in  Figiire  2b, 
the  cold  front  is  that  part  of  the  whole  front  at  which  cold 
air  is  taking  the  place  of  warm  air  at  the  ground,  and  the 
warm  front  is  the  portion  of  the  whole  front  where  cold  air 
is  retreating  with  warm  air  taking  its  place  at  the  ground. 
Since  such  cold  and  warm  fronts  are  accompanied  by  winds, 
and  by  precipitation  in  various  forms --in  fact,  by  all  of 
the  ingredients  that  go  to  make  up  what  one  calls  weather-- 
it  follows  that  the  weather  at  a  given  locality  in  the  mid- 
dle latitudes  is  largely  conditioned  by  the  passage  of  such 
frontal  disturbances.   Cold  fronts  and  warm  fronts  behave 
differently  in  many  ways.   For  example,  the  cold  front  in 
general  moves  faster  than  the  warm  front  and  steepens  rela- 
tive to  it,  so  that  an  originally  symmetrical  disturbance  or 
wave  gradually  becomes  distorted  in  the  manner  indicated  in 
Figure  2c.   This  process  some  times --though  by  no  means  al- 
ways--continues  until  the  greater  portion  of  the  cold  front 
has  overrun  the  warm  front;  an  occluded  front,  as  indicated 
in  Figure  2d,  is  then  said  to  occur.   The  prime  object  of 
what  follows  is  to  derive  a  theory- -or  perhaps  better,  to 
invent  a  simplified  dynamical  model--capable  of  dealing  with 
fluid  motions  of  this  type  that  is  not  on  the  one  hand  so 
crude  as  to  fail  to  yield  at  least  roughly  the  observed  mo- 
tions, and  on  the  other  hand  is  not  impossibly  difficult  to 
use  for  the  purpose  of  mathematical  discussion  and  numerical 
calculation. 

We  begin  with  the  classical  hydrodynamical  equations. 
The  equations  of  motion  in  the  so-called  Eulerian  form  are 
taken : 


Du     3  p    „ 
(3.1,  ^pgv_|£,^,^^^- 
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with  Trr    (the  particle  derivative)  defined  by  the  operator 
+  u  -r —  +  V  -r—  +  w  -:r —  .   In  these  equations  u,  v,  w 


d  t  5x      5y      d  z 

are  the  velocity  components  relative  to  our  rotating  coordi- 
nate system,  p  is  the  pressure,  p  the  density,  PF/  \  etc. 
the  components  of  the  Coriolis  force  due  to  the  rotation  of 
the  coordinate  system,  and  pg  is  the  force  of  gravity  (as- 
sumed to  be  constant).   These  equations  hold  in  both  the 
warm  air  and  the  cold  air,  but  it  is  preferable  to  distin- 
guish the  dependent  quantities  in  the  two  different  layers; 
this  is  done  here  throughout  by  writing  u' ,  v',  w'   for  the 
velocity  components  in  the  warm  air  and  similarly  for  the 
other  dependent  quantities. 

We  now  introduce  an  assumption  which  is  commonly  made 
in  dynamic  meteorology  in  discussing  large-scale  motions  of 
the  atmosphere,  i.e.  that  the  air  is  incompressible.   In 
spite  of  the  fact  that  such  an  assumption  rules  out  thermo- 
dynamic processes,  it  does  seem  rather  reasonable  since  the 
pressure  gradients  which  operate  to  create  the  flov;s  of  in- 
terest to  us  are  quite  small  and,  what  is  perhaps  the  deci- 
sive point,  the  propagation  speed  of  the  disturbances  to  be 
studied  is  very  small  compared  with  the  speed  of  soiand  in 
air  (i.e.  with  disturbances  governed  by  compressibility  ef- 
fects).  It  would  be  possible  to  consider  the  atmosphere, 
though  incompressible,  to  be  of  variable  density.   However, 
for  the  purpose  of  obtaining  as  simple  a  dynamical  model  as 
possible  it  seems  reasonable  to  begin  with  an  atmosphere 
having  a  constant  density  in  each  of  the  two  layers.   As  a 
consequence  of  these  assximptions  we  have  the  following  equa- 
tion of  continuity: 

(2.2)  u  +  V  +  w  =0. 

X    y    z 

The  equations  (2.1)  and  (2.2)  together  with  the  conditions 
of  continuity  of  the  pressure  and  of  the  normal  velocity 
components  on  the  discontinuity  surface,  the  condition  w  =  0 
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at  the  ground,  appropriate  initial  conditions,  etc.  doubt- 
lessly yield  a  mathematical  problem--call  it  Problem  I--the 
solution  of  which  would  furnish  reasonably  good  approxima- 
tions to  the  observed  phenomena.   Unfortunately,  such  a 
problem  is  still  so  difficult  as  to  be  far  beyond  the  scope 
of  known  methods  of  analysis--including  analysis  by  numeri- 
cal computation.   Thus  still  further  simplifications  are  in 
order. 

One  of  the  best-f o'onded  empirical  laws  in  dynamic  mete- 
orology is  the  hydrostatic  pressure  law,  which  states  that 
the  pressure  at  any  point  in  the  atmosphere  is  very  closely 
equal  to  the  static  weight  of  the  column  of  air  vertically 
above  it.   This  is  equivalent  to  saying  that  the  accelera- 
tion terms  in  the  third  equation  of  (2.1)  can  be  ignored, 
with  the  result 

(£.3)  -|£  .  .pg. 

This  is  also  the  basis  of  the  so-called  long-wave  or  shallow 
water  theory  of  surface  gravity  waves,  and  the  author  was 
guided  in  much  that  follows  by  experience  gained  (cf.  [6]) 
in  working  with  such  gravity  wave  problems.   Since  the  ver- 
tical component  of  the  acceleration  of  the  particles  is  ig- 
nored, it  follows  on  purely  kinematical  grounds  that  the 
horizontal  components  of  the  velocity  will  be  independent  of 
the  vertical  coordinate  z  for  all  time  if  that  was  the  case 
at  the  initial  Instant  t  =  0.   Since  we  do  in  fact  assume  an 
initial  motion  with  such  a  property,  it  follows  that  we  have 

(2.4)       u  =  u(x,y,t),    V  =  v(x,y,t),    w  =  0. 

The  first  two  of  the  equations  of  motion  (2.1)  and  the  equa- 
tion of  continuity  (2.2)  therefore  reduce  to 
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u.     +   uu      +  vu      =   -  —   p      +  P/     X 
t  X  y  p  ^x  (x) 

(2.5)  <v.    +UV      +VV     =-tP      +F/V 

\     t  X  y  p   ^y  (y) 

u     +  V      =0, 
X        y 

where  we  use  subscripts  to  denote  partial  derivatives  and 
subscripts  enclosed  in  parentheses  to  indicate  components  of 
a  vector.   The  Coriolis  acceleration  terms  are  now  given  by 


^(x)  ^  2  CO  sin  p 
(2.6) 

P 


•V  =  A 


/    \    =  -2  60  sin  <p 'U.  =   -Au 

when  use  is  again  made  of  the  fact  that  w  =  0.   (The  lati- 
tude angle  Cp,    as  was  indicated  earlier,  is  assumed  to  be 
constant.)   We  observe  once  more  that  all  of  these  relations 
hold  in  both  the  warm  and  cold  layers,  and  we  distinguish 
between  the  two  when  necessary  by  a  prime  on  the  symbols  for 
quantities  in  the  warm  air.   It  is  perhaps  also  worth  men- 
tioning that  the  equations  (2.5)  with  P,  x  and  P,  x  defined 
,    ,  (x)       (y) 

by  (2.6)  are  valid  for  all  orientations  of  the  x,y-axes; 
thus  it  is  not  necessary  to  assume  (as  we  did  above,  for  ex- 
ample) that  the  original  stationary  front  runs  in  the  east- 
west  direction. 

We  have  not  so  far  made  full  use  of  the  hydrostatic 
pressure  law  (2.3).   To  this  end  it  is  useful  to  introduce 
the  vertical  height  h  =  h(x,y,t)  of  the  discontinuity  sur- 
face between  the  warm  and  cold  layers  and  the  height 
h  =  h  (x,y,t)  of  the  warm  layer  itself  (see  Figure  5).   As- 
suming that  the  pressure  p'  is  zero  at  the  top  of  the  warm 
layer  we  find  by  integrating  (2.5): 

(2.7)  ..  p'(x,y,z,t)  =  p'g(h'  -  z) 

for  the  pressure  at  any  point  in  the  warm  air.   In  the  cold 
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Vertical  height  of  the  two  layers 


air  we  have,  in  similar  fashion: 


(2.8) 


p(x,y,zi,t)  =  p'g(h'  -  h)  +  ?^g(h  -  z) 


when  the  condition  of  continuity  of  pressure,  p'  =  p  for 
z  =  h,  is  used.   (The  formula  (2.8)  is  the  starting  point  of 
the  paper  by  Freeman  [3]  which  was  mentioned  in  the  intro- 
duction.)  Insertion  of  (2.8)  in  (2.5)  and  of  (2.7)  in  (2.5)' 
yields  the  following  six  equations  for  the  six  quantities 
u,  V,  h,  u*  ,  v' ,  h'  : 


(2.9) 
(cold  air) 
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u;  +  u'u;  +  v'u^  =  -gh^.  +  /\v' 


(2.10) 


(warm  air)     \ ''t  "^  ^  ^x  *  ^  %  =  '^K   "  ^^ 


X      7^=7 


u;  +  v;  =  0. 


These  equations  together  with  the  kinematic  conditions  ap- 
propriate at  the  surfaces  z  =  h  and  z  =  h' ,  and  initial  con- 
ditions at  t  =  0  would  again  constitute  a  reasonable  mathe- 
matical problem--call  it  Problem  II— which  could  be  used  to 
study  the  dynamics  of  frontal  motions.   The  problem  II  is 
much  simpler  than  the  problem  I  formulated  above  in  that  the 
number  of  dependent  quantities  is  reduced  from  eight  to  six 
and,  probably  still  more  important,  the  number  of  indepen- 
dent variables  is  reduced  from  four  to  three.   These  simpli- 
fications, it  should  be  noted,  come  about  solely  as  a  conse- 
quence of  assuming  the  hydrostatic  pressure  law,  and  since 
meteorologists  have  much  evidence  supporting  the  validity  of 
such  an  assumption,  the  problem  II  should  then  furnish  a 
reasonable  basis  for  discussing  the  problem  of  frontal  mo- 
tions.  Unfortunately,  Problem  II  is  just  about  as  inaccessi- 
ble as  Problem  I  from  the  point  of  view  of  mathematical  and 
numerical  analysis.   Consequently,  we  make  still  further  hy- 
potheses leading  to  a  simpler  theory. 

As  a  preliminary  to  the  formulation  of  Problem  III  we 
write  down  the  kinematic  free  surface  conditions  at  z  =  h 
and  z  =  h'  (the  dynamical  free  surface  conditions,  p  =  0  at 
z  =  h'  and  p  =  p'  at  z  =  h,  have  already  been  used).   These 
conditions  state  simply  that  the  particle  derivatives  of  the 
functions  z  -  h(x,y,t)  and  z  -  h'(x,y,t)  vanish,  since  any 
particle  on  the  surface  z  -  h  =  0  or  the  surface  z  -  h'  =  0 
remains  on  it.  We  have  therefore  the  conditions 
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uh^  +  vh_  +  h.  =  0 

(2.11)  I   u'h  +  v'h  -t-  h   =  0 

X      y    t 

u'h;  .  v'h;  *  h;  =  0, 

in  view  of  the  fact  that  w  vanishes  everywhere.   It  is  con- 
venient to  replace  the  third  equations  (the  continuity  equa- 
tions) in  the  sets  (2.9)  and  (2.10)  by 

(2.12)  (uh)^  +  (vh)y  +  h^  =  0,    and 

(2.13)  [u'(h'  -  h)]   +  [v'(h'  -  h)]^^  +  (h'  -  h),  =  0, 

^  y  t 

which  are  readily  seen  to  hold  because  of  (2.11).   In  fact, 
the  last  two  equations  simply  state  tiie  continuity  condi- 
tions for  a  vertical  column  of  air  extending  (in  the  cold 
air)  from  the  ground  up  to  z  =  h,  and  (in  the  warm  air)  from 
z  =  h  to  z  =  h' . 

We  now  make  the  really  trenchant  assamption,  i.e.  that 
the  motion  of  the  warm  air  layer  is  not  affected  by  the  mo- 
tion of  the  cold  air  layer.   This  assumption  has  a  rather 
reasonable  physical  basis,  as  might  be  argued  in  the  follow- 
ing way:   Imagine  the  stationary  front  to  have  developed  a 
bulge  in  the  y-direction,  say,  as  in  Figure  4a  (see  next 
page).   The  warm  air  can  adjust  itself  to  the  new  condition 
simply  through  a  slight  change  in  its  vertical  component, 
without  any  need  for  a  change  in  u'  and  v',  the  horizontal 
components.   This  is  indicated  in  Figure  4b,  which  is  a  ver- 
tical section  of  the  air  taken  along  the  line  AB  in  Figure 
4a;  in  this  figure  the  cold  layer  is  shown  with  a  quite 
small  height— which  is  what  one  always  assumes.   Since  we 
Ignore  changes  in  the  vertical  velocity  components  in  any 
case,  it  thus  seems  reasonable  to  make  our  assumption  of  un- 
altered flow  conditions  in  the  warm  air.   However,  in  the 
cold  air  one  sees  readily--as  indicated  in  Figure  4c— that 
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quite  large  changes  in  the  components  u,  v  of  the  velocity 
in  the  cold  air  may  be  needed  when  a  frontal  disturbance  is 
created.   Thus  we  assume  from  now  on  that  u' ,  v' ,  h'  have 
for  all  time  the  known  values  they  had  in  the  initial  steady 
state  in  which  v'  =0,  u'  =  const.   The  differential  equa- 
tions for  our  Problem  III  can  now  be  written  as  follows: 

u^  +  uu^  H-  vuy  =  -g[^  h^  +  (1  -  J-)\]    +  Av 

^^•^^^    i      ^t  ■"  ^^x  ■"  ^^y  =  -gfr  ^y  -^  ^^  -  p^^V  ■  ^^ 

(uh)^  +  (vh)   +  h^  =  0, 

in  which  h^  and  h^  are  the  kno\"m  fimctions  given  in  terras  of 
the  initial  state  in  the  warm  air.   The  initial  state,  in 
which  v'  =  V  =  0,  u  =  const.,  u'  =  const.,  must  satisfy  the 
equations  (2.9)  and  (2.10);  this  leads  at  once  to  the  condi- 
tions 


(2.15) 


^„  =  -  —  U  , 

y    g   ' 


h'  =  0 

X 


(h  =-  A, (,.£!),;.  A(£:..  -u). 


V--T^^"  -p^^y  =T^P  "  -^^'   \  =  ° 


for  the  slopes  of  the  free  surfaces  initially.  The  slope  h 
of  the  discontinuity  surface  between  the  two  layers  is  ^ 
nearly  proportional  to  the  velocity  difference  u'  -  u  since 
P'/P  differs  only  slightly  from  unity,  and  it  is  quite  small 
because  of  the  factor  A,  which  is  a  fraction  of  the  angular 
velocity  of  the  earth.  The  differential  equations  for  Prob- 
lem III  can,  finally,  be  expressed  in  the  form: 

f   U^  +  UU^  +  VUy  +  g(l  -  £l)h^  =  Av 

(2-16)         )  o»        .   nt 

Problem  III   \  ^t  "^  ^^x  "^  ^^y  "^  S^^  "  p  ^^y  -  A  (£-  u'  -  u) 

h^  +  (uh)^  +  (vh)y  =  0, 
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by  using  the  formulas  for  h'  and  h'  given  in  (2.15).   We 
note  that  the  Influence  of  the  warm  air  expresses  itself 
through  its  density  p'  and  its  velocity  u' .   The  three  equa- 
tions (2.16)  presumably  have  uniquely  determined  solutions 
once  the  values  of  u,  v,  and  h  are  given  at  the  initial  in- 
stant t  =  0,  and  such  solutions  might  reasonably  be  expected 
to  furnish  an  approximate  description  of  the  dynamics  of 
frontal  motions.  "    Unfortunately,  these  equations  are 
still  quite  complicated,  and  they  cannot  be  integrated  nu- 
merically even  with  the  aid  of  the  most  modern  high-speed 
digital  computers — mostly  because  there  are  three  indepen- 
dent variables. 

Consequently,  one  casts  about  for  still  other  possibil- 
ities, either  of  specialization  or  simplification,  which 
might  yield  a  manageable  theory.   One  possibility  of  spe- 
cialization has  already  been  mentioned:   if  one  assumes  no 
Coriolis  force  and  also  assumes  that  the  motion  is  indepen- 
dent of  the  y-coordinate,  one  obtains  the  pair  of  equations 


(2.17) 


^t  ■"  ^^x  -^  g^^  -  fK  =  ° 


I 


h^  +  (uh)^  =  0 


which  are  identical  with  the  equations  of  the  one-dimensional 
shallow  water  gravity  wave  theory.   These  equations  contain 
in  them  the  possibility  of  the  development  of  discontinuous 
solutions — called  bores — and  this  fact  lies  at  the  basis  of 
the  discussions  by  Freeman  [3,4]  and  Abdullah  [l].   In  such 
one -dimensional  treatments,  it  is  clear  that  it  is  in  prin- 
ciple not  possible  to  deal  with  the  bulges  on  fronts  and 
their  deformation  in  time  and  space,  since  such  problems  de- 
pend essentially  on  both  space  variables  x  and  y.   Another 
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These  equations  are  in  fact  quite  similar  to  the  equa- 
tions for  two-dimensional  unsteady  motion  of  a  compressi- 
ble fluid  with  h  playing  the  role  of  the  density  of  the 
fluid. 
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possibility  would  be  a  linearization  of  the  differential 
equations  (2.16)  based  on  assuming  small  perturbations  of 
the  frontal  siirface  and  of  the  velocities  from  the  initial 
uniform  state.   This  procedure  might  be  of  some  interest, 
since  such  a  formulation  v^ould  take  care  of  the  boundary 
condition  at  the  ground,  while  the  existing  linear  treat- 
ments of  this  problem  do  not.   However,  our  interest  here  is 
in  a  nonlinear  treatment  v/hich  permits  of  large  displace- 
ments of  the  fronts.   One  such  possibility,  to  be  discussed 
in  the  report  by  G.  IVhitham,  involves  essentially  the  inte- 
gration of  the  first  and  third  equations  for  u  and  h  as 
functions  of  x  and  t,  regarding  y  as  a  parameter  and  deriva- 
tives with  respect  to  y  as  negligible  and  assuming  initial 
values  for  v;  this  is  feasible  by  the  method  of  characteris- 
tics.  Afterwards,  v  is  determined  by  integrating  the  second 
equation  considering  u  and  h  as  known,  and  this  can  be  done 
because  the  equation  is  a  linear  first  order  equation  under 
these  conditions.   As  stated  earlier,  this  procedure  fur- 
nishes good  qualitative  results,  but  it  cannot  be  carried 
out  in  all  detail  numerically  because  of  the  fact  that  at  a 
cold  front  one  has,  in  effect,  a  bore  propagating  into  a  re- 
gion containing  no  fluid  and  such  a  situation  cannot  be 
easily  handled.   In  the  following  section  a  different  ap- 
proach to  the  problem  of  approximating  the  solution  of  the 
equations  (2.16)  is  proposed  which  can  be  carried  out  numer- 
ically, but  it  has  the  disadvantage  that  still  more  physical 
assumptions  are  made  and  consequently  the  variety  of  motions 
that  can  be  approximated  is  still  further  restricted. 

§®£ji_2-   Problem  IV:   An  approximate  treatment 
of  the  basic  problem  III 

The  formulation  of  Problem  IV  was  motivated  by  the  fol- 
lowing considerations.   If  one  looks  at  a  sequence  of  weather 
maps  and  thinks  of  the  wave  motion  In  our  thin  wedge  of  cold 
air,  the  resemblance  to  the  motion  of  waves  in  water  which 
deform  into  breakers  (especially  in  the  case  of  waves  which 
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develop  into  occluded  fronts)  is  very  strong.   The  great 
difference  is  that  the  v;ave  motion  in  water  takes  place  in 
the  vertical  plane  while  the  wave  motion  in  air  takes  place 
essentially  in  the  horizontal  plane.   When  the  hydrostatic 
pressure  assumption  is  made  in  the  case  of  water  waves  the 
result  is  a  theory  in  exact  analogy  to  gas  dynamics,  and 
thus  wave  motions  with  an  appropriate  "sound  speed"  become 
possible  even  though  the  fluid  is  incorapressible--the  free 
surface  permits  the  introduction  of  the  depth  of  the  water 
as  a  dependent  quantity,  this  quantity  plays  the  role  of  the 
density  in  gas  dynamics,  and  thus  a  dynamical  model  in  the 
form  of  a  compressible  fluid  is  obtained.   Such  problems 
have  been  much  studied,  and  a  good  deal  can  be  said  about 
them,  principally  because  the  differential  equations  are 
hyperbolic  and  have  only  two  independent  variables.   It 
would  seem  therefore  reasonable  to  try  to  invent  a  similar 
theory  for  frontal  motions  in  the  form  of  a  long-wave  theory 
suitable  for  waves  which  move  essentially  in  the  horizontal, 
rather  than  the  vertical,  plane,  and  in  which  the  waves 
propagate  essentially  parallel  to  the  edge  of  the  original 
stationary  front,  i.e.  the  x-axls .   In  this  xvay  one  might 
hope  to  be  rid  of  the  dependence  on  the  variable  y  at  right 
angles  to  the  stationary  front,  thus  reducing  the  indepen- 
dent variables  to  two,  x  and  t;  and  if  one  still  could  ob- 
tain a  hyperbolic  system  of  differential  equations  then  nu- 
merical treatments  by  finite  differences  would  be  feasible. 
This  program  can,  in  fact,  be  carried  out  in  such  a  way  as 
to  yield  a  system  of  four  first  order  nonlinear  differential 
equations  in  two  independent  and  four  dependent  variables 
which  are  hyperbolic. 

Once  having  decided  to  obtain  a  long-wave  theory  for 
the  horizontal  plane,  the  procedure  to  be  followed  can  be 
inferred  to  a  large  extent  from  what  one  does  in  developing 
the  same  type  of  theory  for  gravity  waves  in  water  (cf.  [6]). 
To  begin  with  it  seems  clear  that  the  displacement  y?(x,t) 
of  the  front  itself  in  the  y-direction  should  be  introduced 
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as  one  of  the  dependent  quantities --all  the  more  since  this 
quantity  is  anyv/ay  the  most  obvious  one  on  the  weather  maps. 
To  have  such  a  "shallow  water"  theory  in  the  horizontal 
plane  requires --unfortunately--a  rigid  "bottom"  somewhere 
(which  is,  of  course,  vertical  in  this  case),  and  this  we 
simply  postulate,  i.e.  we  assiome  that  the  y-component  v  of 
the  velocity  vanishes  f^r  all  time  on  a  vertical  plane 
y  =  o  =  const,  parallel  to  the  stationary  front  along  the 
X-axis  (see  Figure  5).   The  velocity  v(x,y,t)  is  then 


'  h(x,t) 


v7(x,t) 


y='- 


Fig.  5 


assumed  to  vary  linearly'"'  in  y,  and  its  value  at  the  front. 

The  analogous  statement  holds  also  in  the  long-wave  theory 
in  water  in  its  simplest  form. 
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y  =  v7(x,t),is  called  v(x,t).   The  intersection  of  the  dis- 
continuity surface  z  =  h(x,y,t)  with  the  plane  y  =  o  is  a 
curve  given  by  z  =  h(x,t),  and  we  assume  that  the  discontin- 
uity surface  is  a  ruled  surface  having  straight  line  genera- 
tors running  from  the  front,  y  =  W(x,t),  to  the  curve 
z  =h(x,t),  and  parallel  to  the  y,z-plane.   Finally,  we  as- 
sume that  u,  the  x-component  of  the  velocity,  depends  on  x 
and  t  only:  '    u  =  u(x,t).   The  effect  of  these  assumptions 
is  to  yield  the  relations 

(.3.1)  h(x,y,t)  =  J  '  '^I'^'ll   •h(x,t), 

S  -  w (x,t) 

c 

(3.2)         v(x,y,t)=-r — ^^— ^T^-—  .  v(x,t), 

5  -  v7  (x,t) 

as  one  readily  sees  .   In  addition,  we  assiime  that  a  particle 
that  is  once  on  the  front  y  -  v7  (x,t)  =  0  alv;ays  remains  on 
It,  so  that  the  relation: 

(3.5)  v(x,t)  =    ^^   +   uV^^ 

must  hold.   The  four  quantities  u(x,t),  v?(x,t),  h(x,t),  and 
v(x,t)  are  our  new  dependent  variables.   Differential  equa- 
tions for  them  will  be  obtained  by  integrating  the  basic 
equations  (2.16)  of  Problem  III  with  respect  to  y  from  y  =  '^ 
to  y  =  6  --which  can  be  done  since  the  dependence  of  u,  v, 
and  h  on  y  is  now  explicitly  given--and  these  three  equa- 
tions together  with  (3.3)  will  yield  the  four  equations  v;e 
want. 

Before  writing  these  equations  down  it  should  be  said 
that  the  most  trenchant  assumption  is  that  concerning  the 
existence  of  the  rigid  boundary  y  =  o  .   One  might  think 
that  as  long  as  the  velocity  component  v  dies  out  with  suffi- 
cient rapidity  in  the  y-direction  such  an  assumption  would 

The  analogous  statement  holds  also  in  the  long-\'>'ave  theory 
in  water  in  its  simplest  form. 
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yield  a  good  approximation,  but  the  facts  in  the  case  of 
water  waves  indicate  this  to  be  not  sufficient  for  the  accu- 
racy of  the  approximation:   with  water  waves  in  very  deep 
water  the  vertical  velocity  (corresponding  to  our  v  here) 
dies  out  very  rapidly  in  the  depth,  but  it  is  nevertheless 
essential  for  a  good  approximation  to  assume  that  the  ratio 
of  the  depth  down  to  a  rigid  bottom  to  the  wave  length  is 
small.   However,  such  a  rigid  vertical  barrier  to  the  winds 
does  exist  in  some  cases  of  interest  to  us  in  the  form  of 
mountain  ranges,  which  are  often  much  higher  than  the  top  of 
the  cold  surface  layer  (i.e.  higher  than  the  curve  z  =  h(x,t) 
in  Figure  5).   In  any  case,  severe  though  this  restriction 
is,  it  still  seems  to  the  author  to  be  worth  while  to  study 
the  motions  which  are  compatible  with  it  since  something 
about  the  dynamics  of  frontal  motions  with  large  deforma- 
tions may  be  learned  in  the  process.   In  particular,  one 
might  learn  something  about  the  kind  of  perturbations  that 
are  necessary  to  initiate  motions  of  the  type  observed,  and 
under  what  circumstances  the  motions  can  be  maintained. 

Another  objection  to  such  a  theory  could  well  be  raised. 
At  the  very  outset  we  assumed  an  incompressible  fluid,  thus 
ruling  out  at  once  all  thermodynamical  processes.   This 
means  that  we  consider  the  mechanical  sources  of  energy  to 
be  the  essential  ones  in  the  large-scale  motions  under 
study.   Such  a  view  has  been  advocated  by  Abdullah  fl],  who 
supports  it  by  calculations  showing  that  the  over-all  me- 
chanical sources  of  energy  are  of  the  right  order  of  magni- 
tude to  accoixnt  for  the  energy  in  cyclones.   In  the  end,  of 
course,  any  mathematical  theory  which  pretends  to  apply  to 
frontal  disturbances  in  the  atmosphere  can  be  justified  only 
by  its  consequences  in  relation  to  the  observed  facts. 

In  carrying  out  the  derivation  of  the  differential 
equations  of  our  theory  according  to  the  plan  outlined 
above,  we  calculate  first  a  number  of  integrals.   The  first 
of  these  arise  from  (-"^.l)  and  (3.2): 
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/   h  dy  =  y-f—  f     (y  -  vy  )  dy  =  I  h(  6  -    V^  ) 

/  "'^=T^/  (S  -y)  dy  =|v(5  -  v^). 

Prom  these  we  derive  by  differentiations  with  respect  to  x 
and  t  another  set  of  relations: 
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h.  dy  =|h^(6  -  v^)  -fhY^, 
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V   dv=iv(5-V7)+iv 


v..  dy  =iv^(5  -  ^^)  +ivY^, 
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/   ^t  ^y  =l^t(^'  -  7^  -^I^Vt' 


/  h,  dy  =|h,(S  -  >^)  -|hv;,. 

7 

(In  deriving  these  relations,  it  is  necessary  to  observe 
that  the  lower  limit  W  is  a  function  of  x  and  t.)   One  ad- 
ditional relation  is  needed,  as  follows: 
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We  now  Integrate  both  sides  of  the  equations  (2.16) 
with  respect  to  y  from  W  to  8  ,  make  use  of  the  above  inte- 
grals, note  that  u  =  u(x,t)  is  independent  of  y,  and  divide 
by  6  -  Y  •   The  result  is  the  equations 

-     .     -     .   "''Vx  ,  -  ^  7       ^  v^  2kh 

v^  +  uv   +  — +  V.  +  — V/       - 


(3.4)    /t        X     S.y^       t     g.y^rt     5_y^      S-7 

-  2A(u  -  ^  u'), 

uh  +  iiu  -  -^ —  Y)     +  h.  -  -~—    \7.     =  0, 
X    X   c  _  v7   /x    t   S  -  V?   't 


S  -y^   YX    t    5_y^ 


n ' 


with  k  a  constant  replacing  the  quantity  g(l  -  •^) •   These 
equations,  together  with  C?.."^),  form  a  system  of  four  par- 
tial differential  equations  for  the  four  functions  u,  VJ  ,    v, 
and  h.   By  analogy  with  gas  dynamics  and  the  nonlinear  shal- 
low water  theory,  it  is  convenient  to  introduce  a  new  depen- 
dent quantity  c  (which  will  turn  out  to  be  a  "sound  speed" 
or  "propagation  speed")  through  the  relation 

(3.5)  c^  =  I  kh  =  I  g(l  -  ^)h. 

The  quantity  c  is  real  if  p'  is  less  than  p,  and  this  holds 
since  the  warm  air  is  lighter  than  the  cold  air.   In  terms 
of  this  new  quantity  the  equations  (3.3)  and  (3.4)  take  the 
form 

2        , 
u^  +  uu^ -.  2cc^  - -^  y^^  =^  Av 

V   +  uv   =  -  -f2_-  -  2  A(u  -  ^  u'  ) 

2c.  +  cu^  +  2uc   =  °^ 
t      X       X 


(3.6) 


5  -V^ 
^t  ■"  ^^x  =  ^- 
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It  Is  now  easy  to  write  the  equations  (."^.6)  in  the  charac- 
teristic form  simply  by  replacing  the  first  and  third  equa- 
tions by  their  sum  and  by  their  difference.   The  result  is: 

/u^Mu^c)u^-.2Ut^(u-^)c^l-5T;7[Vt^(^^°^7x;=l  ^^' 


(3.7)  ( 


)u^-Ku-c)u^.2Jc,Mu-c)c^j+g^|v^^Mu-c)*?^|  =  |  Av 


-       -         -a-c 
V.  +  uv   =  - 


7 

2  A(u  -  -2^  u'), 


t  X  S   .y^J  '         P 

vn^+uv7^=v. 

As  one  sees,  the  equations  are  in  characteristic  form:   The 
characteristic  curves  satisfy  the  differential  equations 

fTr,\  dx     j^      dx  dx 

(3.8)         dt=^''°'    dt=^-°'    dt=^' 

and  each  of  the  equations  (3.7)  contains  only  derivatives  in 

the  direction  of  one  of  these  curves.   The  characteristic 

dx 
curves  defined  bv  -r-r  =   u  are  taken  twice.   Thus  one  sees 

"  dt 

that  the  quantity  c  is  Indeed  entitled  to  be  called  a  propa- 
gation speed,  and  small  disturbances  can  be  expected  to 
propagate  with  this  speed  in  both  directions  relative  to  the 
stream  of  velocity  u.   (In  the  theory  by  ttTiitham,  In  which 
the  motion  in  each  vertical  plane  y  =  const,  is  treated  sep- 
arately, the  propagation  or  sound  speed  of  small  distur- 
bances is  given  by  Vkh  .   The  sound  speed  in  the  theory 
given  here  thus  represents  a  kind  of  average  with  respect  to 
y  of  the  sound  speeds  of  Whithara's  theory.)   Since  the  prop- 
agation speed  depends  on  the  height  of  the  discontinuity 
surface,  it  is  clear  that  the  possibility  of  motions  leading 
to  breaking  Is  Inherent  in  this  theory. 

Once  the  dynamical  equations  have  been  formulated  in 
characteristic  form  It  becomes  possible  to  see  rather  easily 
what  sort  of  subsidiary  initial  and  bo\indary  conditions  are 
reasonable.   In  fact,  there  are  many  possibilities  in  this 
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respect.   One  such  possibility,  for  v/hlch  numerical  calcula- 
tions are  being  made,  is  the  following.   At  time  t  =  0  it  is 
assumed  that  u  =  const.,  v9  =  0,  h  =  const,  (as  in  a  sta- 
tionary front),  but  that  v7.  =  f(x)  over  a  segment  of  the  x- 
axis .   In  other  words,  it  Is  assumed  that  a  transverse  im- 
pulse is  given  to  the  stationary  front  over  a  portion  of  its 
length.   The  subsequent  motion  is  uniquely  determined  and 
can  be  calculated  niamerically .   Another  possibility  is  to 
prescribe  a  stationary  front  at  t  =  0  for  x  >  0,  say,  and 
then  to  give  the  values  of  all  dependent  quantities  "   at 
X  =  0  as  arbitrary  functions  of  the  time.   One  might  visual- 
ize this  case  as  one  in  ;vhich,  for  example,  cold  air  is  be- 
ing either  added  or  withdrawn  at  a  particular  point  (x  =  0 
in  the  present  case).   This  again  yields  a  problem  with  a 
uniquely  determined  solution,  and  various  possibilities  are 
being  explored  numerically. 

It  v>fas  stated  above  that  the  most  objectionable  feature 
of  the  present  theory  is  the  assu;nptlon  of  a  fixed  vertical 
barrier  back  of  the  front.   There  is,  however,  a  different 
way  of  looking  at  the  problem  as  a  whole  which  may  mitigate 
this  restriction  somewhat.   One  might  try  to  consider  the 
motion  of  the  entire  cap  of  cold  air  that  lies  over  the 
polar  region,  using  polar  coordinates  (9,  <p)    (with  9  the 
latitude  angle,  say).   One  might  then  consider  motions  once 
more  that  depend  essentially  only  on  (f   and  t  by  getting  rid 
of  the  dependence  on  6  through  use  of  the  same  type  of  as- 
sumptions (linear  behavior  in  9,  say)  as  above.   Eere  the 
North  Pole  itself  would  take  the  place  of  the  vertical  bar- 
rier (v  =  01).   The  result  would  again  be  a  system  of  non- 
linear equations--thls  time  with  variable  coefficients.   Of 
course,  it  would  be  necessary  to  begin  with  a  stationary 
flow  in  which  the  motion  takes  place  along  the  parallels  of 
latitude . 

■m 

In  the  niomerical  cases  so  far  considered  we  have  had 
|c!  <  |u|  so  that  even  on  the  t-axis  all  four  dependent 
quantities  can  be  prescribed. 
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All  in  all,  the  ideas  presented  here  and  in  the  follow- 
ing paper  of  Whitham  would  seem  to  yield  theories  flexible 
enough  to  permit  a  good  deal  of  freedom  with  regard  to  ini- 
tial and  other  conditions  so  that  one  might  hope  to  gain 
some  insight  into  the  complicated  dynamics  of  frontal  mo- 
tions . 
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Dynamics  of  Meteorolo.^lcal  Fronts 
by  G.  B.  Whitham 


Sec .  1 .   Introduction 

In  this  paper,  an  alternative  treatment  is  given  of  the 
problem  described  by  Professor  Stoker  in  his  report  "Dynami- 
cal Theory  for  Treating  the  Motion  of  Cold  and  Warm  Fronts 
in  the  Atmosphere."   The  description  of  the  phenomena  which 
we  wish  to  explain  and  the  formulation  of  the  basic  mathe- 
matical equations  are  discussed  in  detail  in  that  report  and 
v;ill  not  be  repeated  here.   The  present  method  starts  from 
the  same  basic  equations  but  differs  in  the  way  in  which 
these  equations  (which  are  too  difficult  to  solve  without 
further  simplification)  are  approximated. 

The  propagation  of  disturbances  on  a  frontal  surface 
appears  to  fall  into  tv/o  parts.   First,  the  Initial  distur- 
bance moves  along  the  discontinuity  surface  with  its  ampli- 
tude (i.e.  the  displacement  of  the  surface  normal  to  its  lin- 
disturbed  position)  Increasing;  then,  at  a  more  advanced 
stage,  the  amplitude  remains  roughly  constant  and  the  cold 
front  begins  to  overtake  the  warm  front,  leading  ultimately 
to  occlusion.   One  may  expect  that  some  qualitative  features 
of  the  first  part  of  the  propagation  could  be  obtained  from 
a  linearized  small-perturbation  procedure,  and  this  has  been 
investigated  by  various  authors  (see  Haurwitz  [5],  p.  307); 
however,  the  occlusion  process  and  the  events  leading  up  to 
it  are  certainly  nonlinear  phenomena.   This  later  develop- 
ment is  considered  here,  and  an  approximate  treatment,  which 
seems  to  be  capable  of  explaining  some  of  the  results  of  ob- 
servation, is  presented. 

It  should  be  remarked  at  the  outset  that,  in  view  of 
the  complexity  of  the  problem,  any  theory  which  includes 
some  of  the  main  effects  or  gives  some  qualitative  insight 
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into  the  problem  is  considered  to  be  of  value,  even  though 
it  may  be  lacking  in  other  respects,  and  may  not  apply  in 
every  situation.   The  present  work  should  be  Interpreted  in 
this  light. 

Sec .  2 .   Approximation  to  the  equations  of  motion 

The  equations  of  motion  for  the  components  u  and  v  of 
the  velocity  in  the  cold  air  and  the  height,  h,  of  the 
frontal  surface  are 


(1)  h^    +    (uh)^   +    (vh)y  =  0, 


(2)  u^    +  uu     +  vu      -   2COV   sin    ©    =   -kh   , 

t  X  y  T'  X 

(3)  V.     +  uv  +  vv      +  2C0  sin    W  {u — i)    =  -kh 

t  X  y  r  p  ;; 


(p   and  p,  are  the  densities  of  the  cold  and  warm  air,  and 
^o     '^l 

u  and  u,  are  the  undisturbed  velocities;  CO  is  the  angular 
o       1  r^ 

Pi 
velocity  of  the  earth,  (D    is  the  latitude,  and  k  =  g(l  -  ■5— )  •  ) 

^o 
In  the  undisturbed  state,  u=u,  v=0,  h=  h  ;  hence,  the 

well-knovm   ;iIarQ;ules  formula   is   obtained: 

dh  n  /-  N     •       /^    Pt  t-It 

(4)  CC  =  2  =  g*^  sin    <p  ,    1   1    _        ) 

^  -'  dy  k  ^    p  ^o'' 

In  general,  u  and  u-,  may  be  functions  of  y  but  it  will  of- 
ten be  convenient  to  take  the  simr)le  case  of  constant  u^  and 

o 

u,  and  to  neglect  the  variation  of  (p    with  y,  so  that  OC  is 
constant  (i.e.  the  cold  air  forms  a  wedge  of  angle  Oc  ) . 

Now,  for  the  situption  which  actually  arises,  OC    is 
small  (a  typical  value  being  1/200)  and  the  theory  developed 
here  is  an  attempt  to  utilize  this  fact.   It  is  assumed  that 
in  the  disturbed  motion  the  rates  of  change  of  the  flow 
quantities  in  the  direction  of  y  remain  small  and  of  the 
same  order  as  OC ;  on  the  other  hand,  h  and  u  may  be  finite 
compared  to  OC ,  and  we  are  interested  in  non-small  values  of 
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them  for  otherwise  we  would  be  reduced  to  something  like  a 

linear  theory.   Then,  in  view  of  the  s.-nall  change  of  h  in 

the  y-direction  and  the  smallness  of  the  Coriolis  terjn,  it 

is  expected  that  the  main  motion  will  be  in  the  x-direction, 

i.e.  the  "crossflow,"  v,  will  be  of  smaller  order  than  u. 

Thus,  approximations  are  made  that  y-derivatives  are  of 

smaller  order  than  x-derivatives  and  that  v/u  is  small. 

More  precisely,  supoose  h  /h  =  0(Cx:);  then,  since 

2CaDu  sin  ^  =  O(kOC)  from  (4),  equations  (2)  and  (l^)  show 

that  Dv/Dt  Si  OC  Du/Dt .   Kence,  unless  v  builds  up  over  a 

large  period  of  time  whereas  u  remains  bounded,  v/u  will  be 

0(OC).   Moreover,  consistent  with  h  /h  C^i  0(0c),  we  should 

expect  u  /u  =  0(0C),  v  /v   =  0(OC).   'Vith  these  approxima- 
y  X  y  x 

tions,  the  ratios  vu  /(vu  ),  vv  /(uv  ),  and  (vh)  /(uh)   are 

y    X     y    X  y     X 

O(a^),  and,  from  (4),  SCOv  sin  ^/kh  =  0(^  ^)  =  O(OC^); 

o   X 

therefore,  neglecting  these  small  terms  (which  are  of  the 

second  order  in  small  quantities),  the  equations  of  motion 

become 


(5)  ■  h^  +  uh^  +  hu^  =  0, 

(6)  u.  +  uu^  +  kh.  =  0, 

(7)  v.  +  uv^  =  260  sin  (D  {-~^   -  u)  -  kh  . 

In  the  first  tv/o  equations,  y  now  enters  only  as  a  parameter, 
but  the  solution  depends  strongly  on  y,  as  vdll  be  seen, 
through  the  conditions  (u  — >  u  ,  h  — *■  h  (y)  as  |x|  — >  co  ) 
under  which  the  equations  are  solved.   l^Tien  u  and  h  are 
known,  v  is  found  from  the  linear  equation  (v). 

Of  course,  in  many  cases  the  approximations  made  above 
may  not  be  true;  for  example,  they  would  not  apoly  when  the 
amplitude  of  the  disturbance  is  small  compared  to  its  wave 
length,  and  thus  could  not  be  used  to  study  the  initial 
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formation  of  frontal  disturbances.   Again,  it  may  be  that  v 
does  grow  continually  until  it  is  no  longer  negligible,   riow' 
ever,  the  approximate  eauations  (5),  (6),  (7)  do  exhibit 
some  features  of  the  motion  which  are  expected  to  be  present 
even  when  other  effects  (neglected  by  them)  also  play  an  im- 
portant role  . 

Sec.  ? .   Solution  of  the  approximate  enuations 

Equations  (5)  and  (6)  have  already  been  studied  in 

2 

great  detail,  since,  on  substituting  c   =  kh,  they  are  iden- 
tical with  the  equations  for  the  velocity  u  and  sound  speed 
c  in  the  one-dimensional  unsteady  flow  of  a  gas  with  ratio 
of  specific  heats  equal  to  2.   With  k  =  gp  they  are  the 
equations  for  waves  on  shallow  water  of  density  p  (see  f6l); 
with  the  present  k  they  have  been  used  by  Freeman  [^]    to 
discuss  waves  on  a  discontinuity  surface  when  there  is  no 
dependence  on  y,  i.e.  the  undisturbed  surface  is  at  a  uni- 
form height  above  the  ground  and  the  disturbance  is  a  f'onc- 
tion  of  X  and  t  only.   In  the  problem  we  are  considering,  -we 
have  essentially  a  Freeman  type  problem  to  solve  in  each 
plane  y  =  constant,  with  the  initial  conditions  depending  on 
the  parameter  y.   Now,  we  consider  the  motion  in  any  such 
plane. 

Introducing  c   =  kh,  (5)  and  (6)  become 

(8)  2c.     +  2uc      +   cu      =0, 

t  XX 

(9)  u.     +   uu      +   2cc      =  0, 

or, 

(.±.   +    (u   +   c)^)(u   +  2c) 
at  -         ox 

=    (u^    +   uu      +   2cc    )    +    (2c.     +   2uc      +   cu    ) 

t  X  X        -  t  X  X 


=  0; 
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hence,  u  +  2c  is  "constant"    on  each  positive  characteris- 
tic curve  dx/dt  =  u  +  c,  and  u  -  2c  is  ''constant"  on  each 
neg:ative  characteristic  curve  dx/dt  =  u  -  c.  ""    This  char- 
acteristic fori-n  of  the  equations  already  enables  any  initial 
value  problem  to  be  solved  numerically  (by  the  well-known 
"method  of  characteristics"),  but  many  important  problems 
can  be  solved  without  numerical  work  by  use  of  the  so-called 
"simple  wave"  solutions.   The  latter  are  solutions  for  which 
either  u  -  2c  or  u  +  2c  is  a  "constant,"  not  only  on  each  of 
the  corresponding  set  of  characteristics  separately,  but 
throughout  the  disturbed  region.   For  example,  suppose  that 
u  -  2c  is  "constant"  in  the  whole  region.   Then,  since  in 
addition  u  +  2c  is  a  "constant''  on  each  positive  character- 
istic dx/dt  =  u  +  c,  it  follows  that  u  and  c  must  be  indivi- 
dually "constant"  on  each  of  these  characteristics.   Kence 
in  the  (x,t)-plane  these  positive  characteristics  form  a  set 
of  straight  lines,  the  slopes  being  determined  by  the  values 
of  u  +  c  on  them. 

An  example  of  the  way  in  which  such  simple  waves  arise 
may  be  seen  in  the  follovdng  problem  (Figure  1,  see  next 
page).   Suppose  that  at  t  =  0,  arbitrary  values  of  u  and  c 
are  given  in  the  range  x-,  ^  x  <  x^,  and  that  outside  this 

range  u  =  u  ,  c  =  c  ,  where  u  and  c   are  "constant."   For 
^       o       o         o      o 

the  moment  we  treat  the  case  in  which  c  is  always  positive 
and  greater  than  u.   Consider  the  negative  characteristic 
BCE  through  x  =  Xp .   It  is  clear  that  alone;  it,  and  along 
any  negative  characteristic  to  the  right  of  it, 

Henceforth,  quotation  marks  around  the  word  constant 
mean   that  the  quantity  is   independent  of  x  and  t;  the 
value  of  the  "constant"  depends  on  y. 

The  case  of  h  =  0,  h  4=  0  at  a  certain  point,  v/ill  arise 
later,  and  it  should  be  noted  that  these  results  do  not 
reduce,  at  c  =  0,  to  u  =  constant  on  dx/dt  =  u,  since  al- 
though c  =  0,  c   =00  and  the  derivation  breaks  down;  in 
fact  the  change  of  u  on  dx/dt  =  u  is  u.  +  uu  =  -kh  . 
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u  -  2c  =  u^  -  2c^  since  all  those  negative  characteristics 


Fig.  1 


meet  t  =  0  at  points  where  u  =  u^,  c  =  c  .   Similarly,  to 
the  left  of  the  positive  characteristic  ACP,  u  +  2c  =  u  +  2c 
In  particular,  then,  for  points  in  EOF  both  the  above  condi- 
tions are  satisfied,  and  we  have  u  =  u  ,  c  =  c   there.   Of 

o       o 
course,  the  region  x  >  Xg  is  undisturbed  until    the  first 
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wavelet,  the  positive  characteristic  BC' :   x  =  x_  +  (u  +  c  )t, 
arrives.   Hence,  the  propagation  of  the  disturbance  to  the 
right  is  confined  between  characteristics  BG  and  CF; 
u  -  2c  =  u   -  2c   in  this  region  so  that  it  is  a  simple  wave 
with  the  positive  characteristics  forming  a  set  of  straight 
lines,  and  u  and  c  "constant"  along  each  of  them.   Similarly 
the  propagation  to  the  left  is  a  si*nple  wave  between  AD  and 
CE,  with  u  +  2c  =  u  +  2c  .   In  the  "interaction"  region, 

O        O  tr>      > 

ABC,  it  is  necessary  to  use  a  rather  devious  numerical 
method.   Two  points  1  and  2  (Figure  2)  are  taken  close  to- 
gether on  the  X-axis  where  u  and  c  are  knovm,  and  short 


->  X 


Fig.  2 

approximate  characteristic  segments  dx/dt  =  u,  +  c,  , 
dx/dt  =  Up  -  Cp,  are  constructed  through  the  points  1  and  2 
respectively  (here  subscripts  denote  values  at  the  corre- 
sponding points);  in  this  way,  a  point  3,  the  intersection 
of  the  two  characteristics,  is  obtained.   The  values  of  u„ 
and  c,  satisfy  the  two  equations  u,  +  2c,  =  u„  +  2c,, 
Up  -  2cp  =  u,  -  2c,  which  are  the  characteristic  relations 
betv/een  u  and  c  on  the  positive  and  negative  characteristics 
respectively;  hence  u,  and  c,  are  easily  found.   The  values 
at  5  can  be  obtained  from  2  and  4,  7  from  4  and  6  and  so  on. 
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The  solution  In  the  region  A.3n  can  then  be  obtained  by  repe- 
titions of  this  process.   When  the  values  of  u  and  c  are 
known  on  BC,  the  simple  v/ave  CPGB  is  obtained  by  drawing 
straight  lines  through  appropriate  points  on  BC  with  slopes 
equal  to  the  values  of  u  +  c  at  these  points.   Similarly  the 
simple  wave  ADEC  may  be  constructed.   When  the  straight 
characteristics  of  the  simple  wave  are  converging,  the  con- 
tinuous solution  is  valid  onlj''  in  the  region  before  they 
overlap.   After  they  overlap  a  discontinuous  jump  of  u  and  c 
across  a  certain  curve  must  be  introduced,  analogous  to  the 
shocks  of  gas  dynamics  and  the  bores  of  shallow  water  theory; 
this  topic  is  postponed  until  Section  4. 

The  above  argument  showing  the  separation  of  an  arbi- 
trary initial  distiorbance  into  two  simple  waves  applies  when 
c  is  sufficiently  large  (i.e.  y  sufficiently  large),  but 
does  not  apply  when  the  plane  y  =  constant  includes  points 
where  the  frontal  surface  meets  the  ground.   Consider  the 
case  when  c  =  0  at  A.  and  C  (Figure  7)    and  suppose  that  there 
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Is  no  fluid  between  A  and  C .   The  locus  In  the  (x,t) -plane 
of  the  point  where  c  =  0  is  dx/dt  =  u  since  equation  (9) 
shows  that,  when  c  =  0,  c^  +  uc^  =  0,  i.e.  c  is  constant  on 
dx/dt  =  u.   Thus  both  characteristic  directions  dx/dt  =  u  +  c 
are  tangent  to  the  locus  of  c  =  0.   Equation  (9)  shows  that 
the  rate  of  change  of  u  along  the  curve  is 

(1^)  \  -^  uu^  =  -kh^. 

Therefore,  if  h^  i  0,  the  characteristic  relations  (u  +  2c  =i 
"constant")  which  degenerate  into  u  =  "constant"  on  c  =  0, 
are  not  satisfied,  and  the  curve  d.:/dt  =  u  is  not  a  charac- 
teristic but  is  an  envelope  of  both  sets  of  characteristics; 
if  h^  -  0,  then  u  remains  constant,  and  dx/dt  =  u  becomes  a 
characteristic.   Again  we  can  argue  that  to  the  right  of  the 
negative  characteristic  through  B,  u  -  2c  =  u   -  2c  ,  and 
there  is  a  simple  wave  on  the  positive  characLristlcs ;  how- 
ever, between  AA'  and  BB',  the  disturbance  would  have  to  be 
determined  namerically. 

One  general  remark  may  be  made  here.   Along  AA',  h  >  0 
and  therefore  from  (10),  u  decreases;  on  CC ' ,  h  <  0,  and" 
therefore  u  increases.   Of  course  the  values  of\  will  be- 
come constant  after  a  certain  time,  but  the  result  shows  at 
least  some  tendency  for  the  gap  to  close  up.   This  is  of 
some  interest  since  it  can  play  a  role  in  explaining  why  oc- 
clusion occurs.   (Another  reason  for  the  overtaking  of  a 
warm  front  by  a  cold  front  will  be  given  in  Section  4.) 

Now  although  it  is  seen  that  an  initial  disturbance  (in 
which  u  and  c  are  both  given  arbitrarily)  does  not  separate 
completely  into  two  simple  waves,  the  author  still  feels 
that  it  is  relevant  to  consider  the  case  of  a  disturbed  re- 
gion which  is  entirely  a  simple  wave  (for  all  y,  not  only 
away  from  the  intersection  with  the  ground).   After  all, 
setting  up  an  initial  value  problem  is  not  necessarily  the 
correct  approach  for  our  problem.   The  frontal  disturbances 
are  initiated  in  some  complicated  way  which  certainly 
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Involves  influences  (the  interference  of  a  mountain  range, 
for  example)  which  are  not  included  in  oiar  ba^ic  formulation 
at  all.   These  may  even  be  more  like  the  piston  problems, 
etc.  of  gas  dynamics  which  are  knowTi  to  produce  simple 
waves.   Moreover,  the  present  theory  makes  no  claim  to  de- 
scribe the  whole  development,  and  is  explicitly  directed  to 
the  later  stages  which  lead  to  occlusion.   We  are  dealing 
here  with  a  fully  developed  disturbance  which  may  have  al- 
ready separated  into  the  two  disturbances  propagating  in  od- 
posite  directions  (a  typical  feature  of  problems  of  wave 
propagation).   In  any  case,  our  object  is  to  see  v/hether  the 
observed  phenomena  are  possible  v/ithin  the  frame v/ork  of  our 
approximate  theory;  for  this  purpose  it  is  natural  to  con- 
sider first  the  simple  waves  which  are  knov/n  to  have  great 
significance  in  other  problems. 

In  particular  the  simple  wave  describing  propagation  in 
the  direction  of  increasing  x  is  considered,  since  fronts 
are  nearly  always  observed  to  be  moving  eastv;ard  (in  the  x- 
directlon)  over  the  United  States.   The  solution  of  (5)  and 
(6)  which  is  chosen,  then,  has  u  -  2c  =  "constant,"  i.e. 

(11)  u  -  2c  =  A(y), 

say,  throughout  the  region.   In  addition,  as  explained  before,  we 
know  that  u+ 2c  is  a  function  only  of  y  on  each  oositive  char- 
acteristic dx/dt  =  u  +  c;  hence,  u  and  c  are  individually 
functions  only  of  y  on  each  of  these  characteristics. 
Therefore,  the  characteristic  equation  may  be  integrated  to 
yield 

(12)  X  =  ^  +  (u  +  c)t, 

where  ^    is  the  value  of  x  at  t  =0.   Now,  the  values  of  u 
and  c  on  the  characteristic  (12)  are  exactly  the  same  as  the 
values  (for  the  same  y)  at  the  point  t  =  0,  x  =  ^  ;  there- 
fore, if  we  suppose,  for  example,  that  c  is  a  given  fvmction 
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C(x,y)  at  t  =  0,  the  value  of  c  on  (12)  is  C(^,y)  and  the 
value  of  u  is,  from  (11),  A(y)  +  2C(^,y).  Thus,  the  com- 
plete solution  is 

c  =  C(5,y), 

(13)  {    ^  =   ^^(y)  "^  2C(^  ,y), 

X  =  "^  +  I  A(y)  +  3C(^  ,y)  I  t. 

(Although  the  arbitrary  function  occurring  in  a  simple  wave 
may  be  specified  in  other  ways,  it  is  convenient  for  our 
purposes  to  give  the  value  of  h,  and  hence  c,  at  t  =  0.)   At 
large  distances  from  the  disturbance  (i.e.  as  ix!  — ^  oo ) ,  u 
and  c  must  have  their  undisturbed  values  '^-.(j)    and  c  (y), 
respectively;  hence,  A(y)  is  determined  by  (11)  to  be 

(14)  A(y)  =  u^(y)  -  2c^(y). 

It  should  be  stressed  that,  in  using  the  solution  (13),  it 
is  much  easier  to  fix  attention  on  individual  planes  y  =  con- 
stant, study  the  solutions  as  functions  of  x  and  t  in  each 
plane  with  y  as  a  parameter,  and  combine  them  finally  to 
give  the  full  picture. 

Turning  now  to  the  equation  (7)  for  v,  it  is  convenient 
(since  u  and  kh  =  c   are  functions  of  ^ ,  y)  to  change  vari- 
ables from  (x,y,t)  to  (^,y,t).   Writing  the  derivatives  in 
the  old  system  on  the  left  and  those  in  the  new  on  the 
right,  we  have 
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where  C,  (x,y)  and  C^{x,j)    denote  3C(x,y)/5x  and  c>C  i.x,j)/  dj, 
respectively.   Equation  (7)  becomes 

l.rj\        <^v       C ^v   _    ___2C___/p     .'^f,  ^ 

^^^^    "aT  -  T-r^c^  -^  -  '  1  +  5c^t(^2  -  ^-  ^^1^ 

+  2C0  sin  ^(-3-^  -  A  -  2C). 
o 

The  details  of  the  derivation  of  the  solution  of  this  equa- 
tion are  not  of  sufficient  importance  to  be  included;  it  may 
easily  be  verified  that  the  general  solution  is 

(18)      v   =  2  J    Q^['^,j)    d^    +    \  C(A'    -   600    sin    </9  ) 


+    (A'    -   200  sin    ^)  5    +  F  j  C^t   +/c^   d^  ,y  j, 

where  P  is  an  arbitrary  function  to  be  determined  from  some 
initial  condition.   For  example,  if  it  is  assumed  that 
V  =  V(x,y)  at  t  =  0  then  F  is  specified  by 

(19)   FJJc^(x,y)  dx,yi  =V(x,y)  -  2  Jc2(x,y)  dx 

-  (A'  -  200  sin  (p  )x. 

The  simple  vi'ave  on  the  other  set  of  characteristics, 
with  u  +  2c  a  function  only  of  y,  and  the  corresponding 
value  of  V  may  be  described  similarly. 

Sec.  4.   The  propagation  of  warm  and  cold  fronts 

The  behavior  of  warm  and  cold  fronts  is  now  discussed 
on  the  basis  of  the  simple  wave  solution  established  in  Sec- 
tion ?.   Across  the  United  States,  the  intersection  of  the 
frontal  surface  with  the  ground  lies  roughly  east-west,  in 
general;  u-,  is  positive  and  u  may  be  of  either  sign  (if 
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u  >  0,  it  must  be  less  than  P-i  i^-i /Pq  ^^   order  that  OC  >  0  in 
(4)).   A  disturbance  on  the  frontal  surface  appears  fre- 
quently to  be  a  northward  push  of  warm  air  and  nearly  always 
is  observed  to  prooaerate  eastward  (increasing  x)  so  that  a 
warm  front  is  followed  by  a  cold  front  as  shown  in  Figure  4 
(see  next  page).   This  being  the  case,  the  siniple  wave  solu- 
tion on  the  "positive"  characteristics  would  appear  to  be 
the  relevant  one,  and  is  the  one  which  will  be  described. 
It  may  be  remarked,  however,  that  as  far  as  the  theory  devel- 
oped above  is  concerned  the  disturbance  could  be  propagated 
either  on  positive  or  negative  characteristics,  or  both. 
The  propagation  velocities  u  +  c  are  positive  and  negative, 
respectively,  for  sufficiently  large  y  (and  hence  large  c), 
but  near  the  intersection  with  the  ground,  v/here  c  is  snail, 
the  direction  of  propagation  would  be  the  same  as  u  in  both, 
cases.   Hence,  if  u  is  usually  positive  (west  v-ind)  there 
would  be  a  tendency  for  the  fronts  observed  at  the  ground  to 
propagate  eastward  in  all  cases.   There  is  some  evidence 
that  u  is  more  frequently  positive  than  negative.   Apart 
from  this  small  argument,  the  explanation  for  the  eastward 
movement  of  the  fronts  must  lie  outside  the  scope  of  this 
theory. 

We  suppose,  now,  that  there  is  a  simple  wave  on  a  set 
of  positive  characteristics  with  u  -  2c  =  u   -  2c   through- 
out, as  given  by  (11)  and  (14).   Immediately  an  important 
result  is  obtained:   At  the  intersection  of  the  surface  with 
the  ground  where  c  =  0,  we  have  u  =  u„  =  u   -  2c  .   If  the 
approximate  values  u  =  constant,  OC  =  h' (y)  =  constant, 
c   =  Vkh  =  s/ocky   ,  are  taken, 


(20)  u^  =  u^  -  2  \/0Cky  , 

so  that  u„  decreases  with  y.   (On  the  cold  front,  (20)  is 
gradually  modified  as  explained  below,  but  remains  true  on 
the  warm  front.)   Moreover,  in  the  plane  y  =  constant,  (1.3) 
gives 
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(21) 


X  =  5  +  j  u^  -  2  \/0C  ky  1 1 


when  C(^,y)  =0;  hence,  since  ^  is  the  x-coordinate  of  the 
front  at  t  =  0,  (21)  shows  that  the  front  has  moved  to  the 
right  a  distance  \u^    -   2  sj  QCkj   Vt  at  ti'ne  t.   This  distance 
decreases  with  y;  therefore,  on  the  basis  of  this  alone,  a 
warm  and  cold  front  system  would  be  distorted  for  successive 
times,  as  shown  in  Figure  4. 

Now,  consider  hov/  the  graph  of  h  against  x,  in  a  plane 
y  =  constant,  varies  with  t.   First  take  y  sufficiently 
large  so  that  C(x,y)  +  0,  for  all  x;  then,  initially,  the 
graph  is  as  shovm  in  Figure  5(1).   At  time  t,,  (13)  shows 
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that  the  value  of  c  at  the  point  x  =  x,  is  equal  to  the 

value  of  c  which  v/as  at  the  point  x  =  ^,  at  t  =  0,  where 

^^  =  X,  -  i  A  +  3C  (  ^,  ,y)  V  t,  .   That  is,  the  value 

c  =  C(^,,y)  has  been  displaced  to  the  right  by  an  arnoiont 

jA+  3C(^-,,y)  ft,.   Since  this  quantity  is  greater  for 
greater  values  of  C,  the  graph  of  h  becomes  distorted  as  in 
Figure  5(ii);  the  "negative  region"  (where  h  <  0)  steepens 
whilst  the  "positive  region"  (where  h  >  0)  flattens  out. 

o  X 

The  positive  region  continues  to  smooth  out,  but,  if  the 
steepening  of  the  negative  region  continued  indefinitely, 
there  would  ultimately  be  more  than  one  value  of  C  at  the 
same  point,  and  the  wave  breaks.  '    Clearly  the  latter  oc- 
ciirs  when  the  tangent  at  a  point  of  the  curve  in  Figure  5 
first  becomes  vertical.   From  (15)  and  (16)  we  have 

-r— -—  — — •  c.>-  (  !^  ,y; , 

ax   1  +  3C^(^  ,y)t   S 

hence,  the  tangent  becomes  vertical  at  t  =  (-3C>-(^,y)) 
The  first  breaking  occurs  on  the  characteristic  ^  =  ^   for 
which  |C^(^  >y)|  is  a  maximum,  and  the  time  of  this  first 
breaking  is 

(22)  t   = 


-1 


^■C^i^^>j) 


m 

At  this  time,  the  continuous  solution  breaks  down  (since  c 
and  u  would  cease  to  be  single -valued  functions)  and  a  dis- 
continuous jump  in  height  and  velocity  must  be  included. 
Such  a  discontinuous  "bore"  propagates  faster  than  the  wave- 
lets ahead  of  it  (the  paths  of  the  wavelets  In  the  (x,t)- 
plane  are  the  characteristics)  in  a  manjner  analogous  to  the 


Throughout  this  work,  since  the  theory  bears  a  close 
analogy  to  the  theory  of  waves  on  shallow  water,  the  ter- 
minology of  water  waves  is  used  in  a  rather  obvious  v;ay. 
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shocks  of  gas  dynajnics.   The  laws  of  propagation  of  a  bore 
are  obtained  by  stipulating  that  the  mass  of  the  fluid 
crossing  it  is  conserved,  and  the  momentiom  of  that  fluid  is 
changed  only  by  the  difference  of  the  pressure  forces  acting 
on  either  side  of  the  bore;  a  certain  amount  of  energy  is 
dissipated  in  turbulence.   (The  details  are  given  in  [6].) 
The  motion  of  the  bore  can  then  be  determined  by  a  method 
which  has  been  used  for  shocks  in  problems  of  gas  dynamics 
(see  [8]);  only  a  qualitative  description  is  given  here 
since  we  are  principally  interested  in  events  where  the 
frontal  surface  meets  the  ground.   (However,  as  suggested  by 
Tepper  fv]  and  discussed  by  Freeman  [5],  the  discontinuity 
in  the  height  of  the  frontal  surface  may  produce  a  sudden 
lifting  of  the  air  above,  resulting  in  the  formation  of  a 
squall  line.)   A  bore  travels  faster  than  the  wavelets  ahead 
and  more  slowly  than  the  wavelets  behind.   In  the  present 
case,  all  wavelets  in  the  negative  region  soon  overtake  and 
run  into  the  bore;  then  the  flow  behind  it  is  in  the  uniform 
undisturbed  state.   At  the  same  time,  the  bore  slowly  pene- 
trates into  the  positive  region  (which  is  continually 
spreading  out);  as  t  — >  oo ,  the  jump  in  height  at  the  bore 
tends  to  zero  and  in  the  remaining  positive  region  h  — >  h  . 
The  (x,t) -plane  is  shown  in  Figure  6  (see  next  page). 

For  a  plane  y  =  constant  in  which  C  =  0  does  occur,  the 
steepening  in  the  cold  front  and  smoothing  out  in  the  warm 
front  take  place  in  a  similar  way  (Figure  7,  see  next  page). 
Again,  the  cold  front  breaks,  and,  in  fact,  since  C  is  pro- 
portional to  the  square  root  of  the  height,  C^  is  infinite 

when  C  =  0  so  that  breaking  commences  immedlatelj''  at  the 
ground.   But,  this  time  there  is  an  important  difference: 
the  breaking  cannot  be  described  by  a  bore  since  there  must 
be  a  flu-x  of  mass  through  a  bore;  hence  the  cold  front  moves 


manner  which  is  difficult  to  treat  theoretically.     However 


over  the  ground,  continually  breaking  and  turbulent,  in  a 
man 


An  approximate  treatment  of  this  problem  has  been  devel- 
oped recently  and  will  be  described  In  a  later  note. 
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we  can  make  some  deductions  about  its  progress  as  follows. 

Before  breaking,  the  velocities  in  the  cold  front  region 

ranp-e  from  u  -  2c   at  the  ground  to  u^  in  the  undistxarbed 
°        o     o  o 

part.   In  breaking,  the  later  regions  carrying  the  higher 

velocities  break  over  the  regions  of  lower  velocity;  hence, 

it  is  certain  that  the  resulting  turbulent  front  moves 

faster  than  u  -  2c   but  slov/er  than  u^ .   Thus,  since  the 
o    o  o 

warm  front  moves  over  the  ground  with  speed  u   -  2c  ,  the 
cold  front  begins  to  overtake  the  warm  front  and  the  occlu- 
sion process  commences. 

The  breaking  of  the  cold  front  starts  at  the  ground  and 
gradually  builds  up  as  more  of  the  fluid  behind  is  included. 
It  would  be  useful  to  have  even  a  rough  guide  to  the  rate  at 
which  it  builds  up  and  overtakes  the  warm  front.   After  a 
time  t,  breaking  v/ill  have  occurred  (in  the  sense  that  the 
tangent  to  the  surface  becomes  vertical)  at  all  points  for 


■r-  1 

Cw(t  ,y)  >  -^  t  (compare  (22)).   Certainly  these  parts 


which 

of  the  fluid  will  have  been  fed  into  the  "breaker,"  and  it 
seems  reasonable  to  take  the  velocity  U  of  the  front  some- 
where between  the  two  extremes  u   -  2c   and  the  velocity  on 

t,  =  ^",  i.e.  u  -  2c^  +  2C(^'"',y),  where  "^  "  is  the  value 

of  ^  for  v/hich  C^(^,y)  =  ^  t.   A  very  simple  choice  v;ould 

be  to  take  the  mean, 

(2:^)  U  =  u^  -  2c^  +  C(^''',y), 

which  starts  at  u   -  2c   (as  it  should)  and  increases  to- 

o     o 

wards  u  -  c  .   (Actually  the  higher  velocities  are  carried 

0      0  o  a 

by  deeper  sections  of  fluid,  hence  the  resultant  velocity 
would  be  weighted  in  their  favor;  on  the  other  hand,  energy 
is  continually  being  dissipated  by  turbulence  and  in  over- 
coming friction  at  the  ground.)   Equation  (2?)  is  extremely 
rough  but  it  should  give  the  correct  order  of  magnitude  for 
the  variation  of  U  with  time,  and  also  something  like  the 
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correct  dependence  upon  y.   In  this  way,  the  intersection  of 

the  frontal  surface  with  the  ground  would  be  modified  as 

shovm  by  the  dotted  lines  in  Figure  4  and  ultimately  the  gap 
closes  up. 

Of  course,  by  the  time  the  gap  closes,  this  theory 
ceases  to  describe  the  motion;  nevertheless  there  is  still  a 
further  prediction  which  can  be  made.   At  the  center  of  the 
occlusion  process  (at  least  initially)  the  velocity  is 
""o  "  2°o  —  ^o  -  2  ^^^y  but  near  the  line  of  the  undisturbed 
front  y  =  0,  the  velocity  is  u^ .   Hence  when  the  gap  closes 
there  would  be  a  tendency  to  leave  a  cyclonic  (anticlockwise) 
rotation  around  the  center  of  occlusion. 

The  decrease  of  the  velocity  at  the  fronts  with  y,  the 
smoothly  propagating  warm  front,  the  turbulent  cold  front, 
and  the  residual  cyclonic  circulation  seem  to  be  in  agree- 
ment with  observation. 

Concluding  remarks 

It  is  hoped  to  do  further  work  in  this  direction  which 
will  include:   (i)  Application  of  the  approximate  theory 
(mentioned  in  the  footnote  to  page  45)  of  waves  breaking 
over  the  ground,  to  the  particular  case  arising  in  the 
"breaking"  cold  front  of  Section  4.   This  will  involve  some 
numerical  integration.   (ii)  Evaluation  of  the  velocity  v 
(obtained  in  Section  3)  in  some  typical  case;  again  numeri- 
cal computation  is  necessary. 

If  the  simple  waves,  discussed  in  detail  here,  are 
found  to  be  inadequate  for  an  understanding  of  the  motion, 
it  will  be  necessary  to  solve  the  approximate  equations  (5), 
(6),  and  (7)  numerically.   It  may  be  remarked  that  the  way 
in  which  such  a  solution  would  be  carried  out  is  quite  stan- 
dard and  well-known,  and  should  not  give  rise  to  any  addi- 
tional difficulties;  the  difficulty  of  the  wave  breaking  at 
the  intersection  of  the  discontinuity  surface  with  the 
ground  on  the  cold  front  side  is  still  expected  to  arise 
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